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DROUGHT FLOWS AND RECEIVING WATER ASSESSMENT 



Lloyd A. Logan 



ABSTRACT 

n.^.^.f-""^"^^^^ water quality objectives have been established for the 
protectTon of aquatic life and for recreation The hvdro ooic 
characteristics have don,inant effects on the abi ity of the^ treamflow 
ob1ecti:e ''%rr''"'^"" '° "^^' '"' P-°«^"cial water qua t^ 
cr ter a for ^c^I'"" '"^ f^^^^^^^^ are used to establish design flow 
criteria for receiving water assessments. Further the effects nf 
drought on streamflows are considered by evaluating performances in 
satisfying design flows with regards to possible impacC oHstp '" 
quality protection A statistical .ethod ha beeZ eS to cha acterize 

ow "r^ s Z :''f '"°"' V'"' '"^° probable critic dough' 
con^H»rfn^f=M P^'^f"™""^ °f the rivers were examined fay ^ 
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INTRODUCTION 

Several rivers in Ontario' are used tO' asslmillate treated 
wasteloads. Mian made wastes include point source discharges from 
municipal and industrial plants and diffuse sources in runoffs. These 
waste sources may have miajor impacts to the water quality. Safe water 
quality is necessary for a viable aquatic life and to provide stream 
conditions which are free of risk to human health. 

The Ontario Ministry of the Environment has the responsibility to 
protect the environment from degradation (Government of Ontario 1975, 
1980, 1980). Guidelines and objectives have been established for water 
quality protection (Ontario Ministry of the Environment 1979, 1980, 
1983); particularly, to protect and improve the quality of streams and 
lakes that have been degraded and to minimize the risk of failure in 
flows for wastewater assimilation. 

The hydrological conditions influence the capability of the river 
to assimilate wastes. Droughts, in particular, are critical to waste 
assimilation. Similarly, the physiographic and hydrogeoTogic 
characteristics vary in nature in Ontario (Chapman and Putman 1966) and 
subsequently influence the streamflow yields which are natural sampling 
of the water ava i 1 ab i 1 i ty . 

The provincial streamflow yields belong to a spatial and definite 
probability distribution. In addition, the time series of extreme low 
flows belong to a separate and different probability distribution. A 
number of non-parametric probability distributions have been 
investigated for extreme low flows (Gumbel 1958, 1963; Matalas 1963). 
A comparative evaluation of probability distributions for extreme low 
streamflows in Ontario were investigated by Condie and Nix (1975). The 
Ministry of the Environment has analysed extreme low flow time series 
in Ontario to establish design flows for wastewater assimilation (OMOE 
1973, 1974, 1975, 1976, 1977). Further investigation has evaluated the 
suitability of an acceptable risk of failure to set guidelines for 
design flows (Logan 1982, 1984). 

The objectives of this paper are to evaluate ex:trem:e low' flows 

with respect to meeting the needs of instream quality assessment and to 
evaluate the perforraiance of the river system in terms of the streamflow 
reliability, resiliency and vulnerability. 
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DROUeHT FLOW DESIGN 

The Provincial streamfloiw yield, being spatial in nature,, is 
assumed to be norniany distributed: N [p, cr] with confidence bounds 
about the population average,,|j , defined as: 



[1] X - t S S M S X + t^ S 



where, X and S'x denote the samiple average and standard error of 

average respectivily; t^ denotes the standard deviate with a the 
significance level. 

The extreme low flow event, on the other hand, is defined as the 

iiiniroum consecutive ra-day average low flow in any given year: 



[2] X (ii)j - min [X (ni).^.] 



where, X(m;)j and X (m)i_j denO'te the consecytive m-day 

average extreme low flow in year, j and the consecutive m-day moving 

average. flow at period, i, in year, j, respectively. 

The time series of extreme low flows, X (m) j , j - 1,2, ..,N, at a 
given site, is assumed to be Gumbel Extremal Type 3 distributed 
(Gumbel, 1963): 



[3] f(x) = [^] [^V" exp -[^] 





where, Xq is the .characteristic drought flow; a is the scale 
.parameter, (equivalent to standard deviation); and e is the lower 
boundary value; and with 



X = X (m). and x = X (m). . are defined 

J II ■» J 
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Integration of eq. [3] from e to x produced the cumulative 
probability distribution: 



[4] P(x) ^ e"^ 
.where, y - 



X - e l 



"X ~ 




A drought flow design, mQj, can be specified from the fitted 
probability model, eq. [3J. The risk of failure T (return period) is 
set by the decision-maker or from guidelines. An acceptable design 
condition is subjected to a specified water quality objective (OMOE 

1980i). 

Parameter Uncertainties 

Uncertainties are imbedded in the sample estimates of the 
probability model parameters, XQ,a, and e (eq. [3]). A computer 
intensive non-parametric technique (Efron and Gong 1983; Logan and 
Unny 1985) Is used to investigate these uncertainties. 

The procedure involved fitting the extreme low flow events, 
Xj J = 1, 2, .., N, to the probability roodel , to determine estimates 
of Xq, a, and e. By applying the Monte Carlo resampling procedures 
implicit in the Boothstrap's approach, several realizations of Xq, a, 
and e can be determined. Frequency analyses of the realizations of the 
parameter sets provide the joint empirical distributions which describe 
the uncertainties associated with these parameter estimates. 

Waste! pad 

Site specific design flows for wasteloads have been investigated 
by Biswas and Bell (1984). They have considered a design flow to 
satisfy chronic and acute values of pollutants. For protection against 
chronic value, the wasteload, W^* at a site, is proportional to the 
total flow times the concentration of the pollutant: 



[5] W a (raQ-r + Q ) C 

c - T ^c c 



where, W^ is the chronic wasteload; mQj is the design low 

streamflow; Q^ and C^ are the effluent flow and concentration, 
respectively. 
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The non-conservative pollutants, example, BOD and other wastes, 
depend on pH and temperature for microbial breakdown of the biological 
matters; hence to determine the chronic wasteload, W^ » an empirical 
water quality model analysis is required to relate flow and waste (OMDE 
1976, 1977). For this study it is assumed that once the W^ is 
determined, this effluent disposal can be controlled, in time, so as 
not to be exceeded. There is, however, chances of violation of the 
specified water quality objective due to stochastic variability in 
stfeamflow. This implies that the non-violation of the water quality 
objective is conditional on the ability of the river to satisfy the 
design flow requirement at all times. Logan (1984) investigated risks 
of failure in water qyality considering stochastic variabilities in 
both streamflow and quality. This paper will consider only streamfTow 
performance with implication and inferences on likely risk of failure 
in water qua! ity. 



STREAMFLOW PERFORMANCE 



If the streamflow events over time period of N years are such 
that the targetted drought flow design, mQi-, is satisfied at all 
times and the chronic wasteload, W,^, discharged is not being 

exceeded,, then the streamflow system is reliable.. 

Hashimoto et al (1982) developed concepts for evaluating system 
performance based on reliability, resiliency and vulnerability. 
Fiering (1982) quantified the system resiliency by considering how 
quickly a system can recover from failure once a failure occurred. He 
further expanded this to the severity of failure once failure occyrred 
which he defined as vulnerabn ity. 

Rellabiility of the system is the probabil ity,a, that the system 
is in a satisfactory state, S,: 

.[6] a = P [Xt £ S] 

where, X^^ denotes the streamflow at time t. 

The degree to which the system recovers once a failure occyrred 
is the resil lency: 

[7] ■ Z = P[Xt+;L, ^ S' » '^'t ^ F] 

where, F denotes that the system is in a failure state. 

The system operates in transition from a satisfactory to a failure 
state and vice versa. 

Let Z^ = 1, if Xi^ E S; and Z^ = 0, if X^ s F; hence 

the fraction of the time the system Is satisfactory is: 

[8] Ts = (1/n) Z Zt 

t = 1 
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In the long run as n -+ «©, satisfaction approaches the system 
rel iabil ity; 

n 
[9] Lim 1/n Z Z^ = a 

n -► «o t=l 

Consider the process of transition: 

Let Wt - 1, if Xt E S and Xt+i s F; and W^ = 0, 
otherwise. 

Similarly, in the long run as n -* «», the average value of W^ 
approaches the probabil ity,p , of the system being in the state, S, and 
transferring to the failure state, F, in the following period: 

n 
[10] p = lim 1/n E Wt = P[Xt £ S; Xt+i z F] 
n -+- <» t=l 

where, p is the probability in transition. The average time in 
temporary resident in the unsatisfactory state, F, during a n-period 

is: 

[11] Tp = (total time spent in F)/( number of times F occurred) 
Hence: 



n • n 

[12] Tp- [1/n E (1=Z^)] / [1/n E W^] 

t=l t=l 

^ (1- C£)/p 



By definition, resiliency,!, is the reciprocal of E[Tp]; therefore 

with n -* «•: 

[13] IT = p/(l=a) = P[Xt E S; Xt+i z F]/P[Xt £ F] 

Further, as n -♦ » the expected number of transition from S to F is 
equal to the expected number of transition in the reverse direction 

[14] P[Xt E S; Xt+i e F] = P [Xt t F; Xt+i s S] 

This implies that If is equivalent to the average probability of 
recovery from the failure, F, in a single time step: 

[15] Z = P[Xt E F; Xt+i E S]/P [Xt £ F] 

= P[Xt+i £ S II Xt E F] 
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Vulmerabiil iity is a, measure of the severity of a failure in a 
system and, as such, defines the magnitude of the failure. This could 
either be in termis of intensity and duration. Although a system may 
not likely to fail it is, however, still vulnerable. 

Let X^ be the discrete failure-flow of the system output, that 
is, X-^ E F; and let Df^^ denotes the severity or the degree of 
failure assigned to each failure. In drought flow periods, a low 
intensity long duration failure is considered to be as imiportant, as a 
high intensity short duration failure, given that an equivalent short 
fall quantity Is experienced for each failure occurrence. Therefore, 
the degree of failure, Dj^, due to the occurrence of a failure on any 
single day, t, is expressed in terms of the deficit between the system 
requirement and the prespecified design flow, rapj, desired for the 
wastewater assimilation, that is, using both intensity and duration to 
define failyre gives: 

[16] Dj^ - (mQT - X|^) A t/A 

where, A t is the period between t and (t + 1), A is a constant used 
to standardize ID|^; and k is the time step of the severe failure-floiw. 

Further, define h^ as the probability that a discrete failure, 
X\, is the most severe outcome that occurs within the set of' failure 
states. Then vulnerability of the system, v, is specified by: 

L 
[17] V = I hk D'l, 

k = 1 



where, k = 1, 2, ..., IL is the number of severe fai lure-flow events. 

Trade-offs between system resiliency and vulnerability is an 
important factor in decision making with regards to operational 

pol icies. 



APPLICATION 



There are approximately 505 gauged streamflow locations in 

Ontario. These streams are surveyed by the Water Survey of Canada; 
with published records issued annually (WSC 19XX). The stations are 
identified by Drainage Basin Boundary as shown by examples in Figure 
1. 

Sets of average annual flows at 233 natural and 272 regulated 
stations indicated Normal Distributions as shown in Figure 2(a); with 
the provincial yields and standard errors of 0.0111 + 0.0002 m^/s/km^ 
and 0.0117 ± 0.0002 m^/s/km^ for the natural and regulated sites, 
respectively. 
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Figure 1. Ontarioi: Examples of Streamfllow Station Location, 
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Figure 2. Distribution of Streamflow Yeilds in Ontario. 
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Frequency analysis of extreme low flows at gauged stations with 
records greater than lO' years prodluced results at 76 natural and 111 
regulated stations (OMOE 1973, 1974, 1975, 1976, 1977). Spatial 
frequency analysis for design low flows: 7Q2, 7Qiq and 7Q2o 
demonstrated skewed distributions (Figure 2(b)). The skewness 
increases negatively with an increase in the return period. Further 
spatial analysis of design low flow yields, on a regional basis, showed 
noted difference in average yields between the southern and the 
northern regions of the province (Table 1). The southwest region has 
yields which are significantly lower than the provincial average 
yields, implying greater probability for susceptibility to severe 
drought. On the other hand, the northern, northeast^ and central 
regions appear to be less susceptible to drought; in that, they are of 
expected higher streamflow yields which appeared to be consistent with 
the natural yields of the hydrogeological condition of those areas. 
The spatial yield values are shown only for comparative purposes and 
should not be used in design. 

Characteristic drought flows frequency curves at selected 
streamflow locations for defined consecutive m-day drought sequences, 
calculated for the time period 1960-1984, are shown in Figure 3(a). 
The rivers in the north and central Ontario (example, Kaministiquia and 
North Magnetawan) display higher drought yields as compared to the 
sooth (example, last Humber and South Nation). Other examples of site 
specific drought flow frequency curves are shown in Figures 3 (b). 
Additional frequency curves for drought flows can be derived with the 
use of the values listed in Table 2. Unlike the regional low flows, 
the stations' drought flow frequency curves can be used for design 
purposes. 

The parameters: characteristic drought flow, Xq, and scale 
parameter, a, are examined for uncertainties using the Boothstrap 
computer intensive statistics. In this example, the minimum daily 
flows are fitted to the Gumbel (Fisher-Tippett Type 1) distribution. 
The Monte Carlo resampling approach in 1000 Boothstrap iterations 
produced the joint distributions for Xq and a as shown in Figure 4 
for selected rivers. These joint distributions are bounded and 
appeared to be inherently Normal. Each pair of parameters (Xq and a) 
of the empirical distributions can be used to generate a series of low 
flows which are utilized in further risk analysis of a separate 
study. 

The reliability of a river is its capability to satisfy a target 
or design f 1 ow , mQf , at a site on a daily basis, wh i 1 e res 1 1 1 ency 
measures are the probability that the river will recover from failure 
in time immediately following the time the failure occurred. In fact, 
a failure is said to occur when the output violates a specified 
performance threshold. If a failure in mQT occurred at a control site, 
it is then implied that a water quality violation will occur. If the 
failure is of a considerable duration, in terms of consecutive days of 
non-exeedence in mQj, then the impact on water quality violation 
could be insurmountable. 

The Ministry of the Environment's Policy is a consecutive 7-day 
average Tow flow at a 5 percent risk of failure (7Q2o). It is 
implied that non-exceedences of the design flow will cause water 
quality violation in BOO and other wastes. Vulnerability gives a 
measure of the severity or magnitude of this failure* 



TABLE 1: AVERAGE YIELD FOR DESIGN FLOW: ONTARIO 
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.10-3] 
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0.2 


P.I 


7Qio 
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Figure 3. Drought Flow Frequency Curves 
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TABLE 2: DROUGHT FLOWS AT SPECIFIED RETURN PERIODS 



DROUGHT 
SEQUENCE/ 


KAM. 


PIN. 
[Flow 


H-MAG. 
Yield (b3 


NOTT. 
/S/Krm2).10- 


E-HUM. 

'1 


S-NAT. 


1-Day: 


7Q2 


4.60 


1.44 


2.55 


1.89 


0.65 


0.24 




7Q10 


1.68 


0.94 


0.90 


1.22 


0.25 


0.11 




7Q20 


1.40 


0.84 


0.60 


1,12 


0.14 


0.10 


3-Day: 


7Q2 


4.84 


1.40 


2.62 


1.96 


0.72 


0,27 




7Q10 


2.74 


0.94 


0.96 


1.46 


0.42 


0.11 




7Q20 


2.29 


0.84 


0.63 


1,30 


0.36 


0.10 


7-Oay: 


7Q2 


5.77 


1.40 


2.78 


2.03 


0.77 


0.28 




7Q10 


3.54 


0.94 


1.06 


1 ., 52 


0.47 


0.15 




7020 


3.06 


0.84 


0.69 


1.37 


0.41 


0.14 


15 -Day: 


702 


6.44 


1.41 


2.98 


2.17 


0.84 


0.33 




7Q10 


4,15 


0.95 


1.28 


1.63 


0.52 


0.17 




7Q20 


3.70 


0.85 


0.92 


1.49 


0.45 


0.16 


30-Day: 


7Q2 


7.39 


1.44 


3.12 


2.36 


1.00 


0.42 




7Q10 


4.55 


0.99 


1.54 


1.75 


. 64 


0.18 




7Q20 


3.94 


0.90 


l.,12 


1.62 


0.54 


0.15 
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The characteristic drought flows relate to the greatest impacts a 
defined drought sequence may have at a site should an impending drought 
occur. The performance of a river is evaluated by exainining the 
reliability, a, resiliency, Z, and vulnerability, v. Table 3 shows the 
performance evaluations based on the characteristic drought flows at 
selected rivers in Ontario. Graphical descriptions are shown in Figure 
5. The vulnerabilities appeared to be marginal for defined drought 
sequences up to the 30-day average. These are due to the Tow flow 
deficits experience during failure days. The resiliencies are, however, 
very critical, resulting in extended duration of consecutive days of 
failure, once failure occurred. Eventually, these failures will have 
a major impact on the water quality. The stations: Nottawasaga, 
Humber and South Nation rivers displayed similar performances, implying, 
high reliability, marginal vulnerability and concerned cases of 
resiliency. The Pinelmuta and Kaministiquia rivers, on the other hand, 
demonstrated a greater measure.=of vulnerability, while the resiliencies 
remain critical. The station North Magnetawan showed high 
reliabilities and with resiliency values ranging from 0.027 to 0.026; 
implying a range of 37 to 39 consecutive days of failure, for defined 
l-day to 30-day average drought, once a failure does occur. The Nith 
River showed a similar performance of implied duration of 25 to 36 
consecutive day of failure once a failure occurred; and of lesser 
reliability. 

For design purposes in satisfying wastewater assessments, the 
magnitude of acceptable performance of the river, mainly resiliency, 
should be enhanced. Performance evaluations based on drought flow 
designs: 7Q2, 7Qio and 7Q20 demonstrated significant 
Improvement on satisfying performance requirement. Table 4 shows 
values for ot, If and v at selected sites. Figure 6 displays the 
characteristic performance curves. Reliabilities in all cases are 
high, resiliencies are of acceptable value, indicating reduction in 
duration of consecutive days of failure once failure occurred. In 
general, resiliency and vulnerability in comparison to the 
characteristic drought flows are at higher efficiency, but are at lower 
improvement between policies 7Q2 and 70io* however, performances 
between policies IQiq and 7Q20 are at acceptable levels 
implying high reliability, moderate resiliency (reduction of 
consecutive days of failure once failure occurred) and negligible 
vulnerability. Suitable management practice could best be evaluated 
for meaningful trade-off between Z and v. 
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TABLE 3: PERFORMANCE 


EVALUATION 


OF RIVER SYSTEM UNDER CHARACTERISTIC DROUGHT FLOWS 


STREAMFLOW 


SEQUENCE 


FLOW f 


RELIABILITY 


RESILIENCY 


VULNERABILITY 


CONSECUTIVE 


STATION 


(-day) 


(mVs/K»2) 
lOr'3 




(T-M 


{m^} 


FAILURE DAYS 


KAMINISTIQUIA 


1 


4.78 


0.98 


0.281 


0.11 


>*. 


(02AS010) 


3 


5.48 


0.95 


0.138 


0.28 


7, 




7 


6.47 


0.91 


0.055 


2.23 


18. 




15 


7.19 


0.88 


0.046 


3.98 


22. 




30 


8.29 


0.84 


0.043 


5.83 


23. 


HORTH 


1 


3.08 


0.93 


0.027 


0.33 


37. 


MAGNETAWAN 


3 


3.15 


0,92 


0.,028 


0.32 


36. 


(02EA005) 


7 


3.30 


0.91 


0.028 


0.35 


36. 




15 


3.49 


0.91 


0.026 


0.42 


39. 




30 


3.64 


0,90 


0.026 


0.45 


39. 


NOTTAWASAOA 


1 


2.02 


0.95 


0.1012 


0.04 


10. 


(02ED003) 


3 


2.07 


0.94 


0.094 


0.05 


11. 




7 


2.15 


0.92 


0.078 


0.05 


13. 




15 


2.31 


0.89 . 


0.056 


0.11 


18. 




30 


2.55 


0.83 


0.034 


0.33 


30. 


EAST 


1 


0.74 


. 94 


0.107 


0.001 


9. 


HUMBER 


3 


0.84 


0.92 


0.075 


0.002 


13. 


(02HC009) 


7 


0.85 


0.91 


0.055 


■ 0.003 


18. 




15 


. 94 


. 87 


0.048 


0.004 


21. 




30 


1.10 


0.80 


0.037 


O.OOB 


27. 


SOUTH 


1 


0.30 


0.94 


0.060 


0.06 


17. 


NATION 


3 


0.31 


0.93 
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Figure 5. Performance Evaluation o^f River System 
Under Characteristic Drought Flow(Xg). 
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TABLE 4: RIVER SYSTEM PERFORMANCE: DESIGN FLOW 
Reliability: a = P [Xt S JQj] 



POLICY 


KAMINISTIQUIA 
(%) 


N-MAGNETAWAN 
(%) 


NOTTAWASAGA ' 
(%) 


7Q2i 


99.7 


100.0 


99.6 


7QiB 


99.6 


99.8 


99.3 


7Q2 


93.9 


94.3 


94.3 



Resiliency: I - P[\+i S 7Qjl!Xt< 7Qx] 



POLICY 


KAMINISTIQUIA 


N-MAGNETAWAN 


NOTTAWASAGA 


7Q,, 


0.56 , 


1.00 


0.21 


7Qii 


0.39 


0.17 


0.20 


7QI2 


0.10 


0.03 


0.10 



Vylniiraibil ity; V = Z h^ U\^ 



POLICY 


KAMIINIISTIQUIA 


N=IMAGNETAWANI 


NOTTAWASAGA 


7Q20 
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0.00 


0.003 


7Qio 
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0,003 


0.003 
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0. 08 


0.040 



s 



UJ 
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Figure 6. Performance Evaluation of River System 
Under Drought Flow ID'esigini(7Qi^) 
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CONCLUSION 

A frequency anaylsis of annual average streaniif low yields for 
the province has identified Normal Distributions for both the 
natural and regulated streamf lows . The nion-parametric probability 
model, Gumbel Extremal Type 3 Diistribution , is adopted for 
evaluating extreme drought flows in Ontario. Characteristic 
drought flows have been derived for gauged stations with 
corresponding frequency curves for specified return periods. 
Probability distribution, on a spatial basis, for the consecutive 
7-day drought flow designs 7Q2 , IQia s^d 7Q20 were shown to be 
negatively skewed, indicating non-normality. 

A computer-intensive statistic, using the Boothstrap 
technique, has been used to develop site specific empirical 
non-paramietric joint distributions for the characteristic drought 
flow, Xq, and the scale parameter, a, of the Gumibel Type 1 Model. 
The distributions are bounded, but adequately described the 
uncertainties associated with these parawiieters. 

Perforfnance evaluations for selected rivers in Ontario, 
based on the characteristic drought flows and assigned drought 
flow designs 72, IQiq and 7Q2i have demonstrated satisfactory 
performance statistics a, I and v. The reliabilities were 
satisfactory in all cases except that reliability by itself is not 
a sufficient statistic to evaluate the river system performance. 
The resiliency and vulnerability, in addition, showed m_arked 
changes with differing policies. 

In the case of the characteristic drought flows, the vulnerability 
increases and the resiliency decreases with an increase in days of 
the defined drought sequence. This imiplied that although the 
magnitude of flow to recover from failure may not be large, the 

resiliency nonetheless is critical; due mainly to the length of 
duration of the consecutive days of failure once a failure 
occurred. The resiliency factor is critical to both the southern 
and northern rivers under characteristic drought flow conditions. 

iMianagement decision was examined utilizing drought flow 
designs projected to meet specific needs of wasteload 
assimilation. Significant improvements in the performance 
evaluation of the rivers were achieved. Less improvement was 
however observed for system performance between the policies JQg 
to 7Q.10, as compared to values between policies 7Qin. and 7Q120 
where a greater degree of tolerance were achieved. That is, the 
degree of vulnerability was reduced markedly and the duration of 
consecutive days of failure once a failure occurred, reduced to 
tolerable durations. Some differences were observed between 
streamf low stations in the north where both vulnerability and 
resiliency would be of concern as compared to the stations in the 
south, where resiliency would have the greater imipact on the 
system especially to water quality. The current Ministry's 
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guiideline, on drought flow design, 7Q2o. demonstrated exceinent 
system perforniance evaluation at most sites. In that the 
streamflow performance for wastewater assimilation, showed high 
reliability and tolerable resilience. Greater susceptibility to 
shock of failures and damages to the river once failure occurred 
are therefore more critical especially in the characteristic 
drought flow ranges. Suitable management decision on design 
policy for wastewater assessment could avert likely potential 
imipacts of drought to water quality. 
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